clear;
M=[1 1 0 0 0 0 0 0 0 0 0 0;
 1 6 1 0 0 0 0 0 0 0 0 0;
 0 1 6 1 0 0 0 0 0 0 0 0;
 0 0 1 6 1 0 0 0 0 0 0 0;
 0 0 0 1 6 1 0 0 0 0 0 0;
 0 0 0 0 1 6 1 0 0 0 0 0;
 0 0 0 0 0 1 6 1 0 0 0 0;
 0 0 0 0 0 0 1 6 1 0 0 0;
 0 0 0 0 0 0 0 1 6 1 0 0;
 0 0 0 0 0 0 0 0 1 6 1 0;
 0 0 0 0 0 0 0 0 0 1 6 1;
 0 0 0 0 0 0 0 0 0 0 1 1;];
B(1,1)=2.0*1/(1+5*5);
for i=1:10
 B(i+1,1) = 8.0*1/(1+(i-5.5)*(i-5.5));  
end
B(12,1)=2.0*1/(1+5*5);
C=inv(M)*B;